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The inverse eigenvalues of the Dirac operator in the Schwinger model satisfy the 

(N 

' same Leutwyler-Smilga sum rules as in the case of QCD with one flavor. In this 

paper we give a microscopic derivation of these sum rules in the sector of arbitrary 
I topological charge. We show that the sum rules can be obtained from the clustering 

(N 
O 

vn 
o 



property of the scalar correlation functions. This argument also holds for other 
theories with a mass gap and broken chiral symmetry such as QCD with one flavor. 
For QCD with several flavors a modified clustering property is derived from the low 
energy chiral Lagrangian. We also derive sum rules for a fixed external gauge field. 
Both types of sum rules are obtained from the bosonized version of the Schwinger 
model. 

^ I In the sector of topological charge u the sum rules are consistent with a shift 



Oh! 



of the Dirac spectrum away from zero by average level spacings. This shift is 
also required to obtain a nonzero chiral condensate in the massless limit. Finally, 
we discuss the Dirac spectrum for a closely related two-dimensional theory with a 
gauge field action that is quadratic in the gauge field and has applications to the 
quantum Hall effect and d-wave super-conductors. 
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I. INTRODUCTION 

The SchwiiiErer model hj, electrodynamics in 1 + 1 dimension is one of the few quan- 

nnn 

tum field theories that can be solved analytically y, 1^- Since the model contains 
nonperturbative features that are also found in QCD, such as a nonzero chiral conden- 
sate and instantons, it has added significantly to our understanding of gauge theories 
in general. However, most interesting quantum field theories cannot be solved analyt- 
ically, and one has to rely on approximate techniques to investigate them. For non- 
perturbative phenomena, two approaches which both have been tested in the Schwinger 
model, have been very successful: lattice simulations and chiral perturbation theory. In 
jarticular, a larg e number of lattice studies have benefited from the Schwinger model 



jarticuiar, a large number oi la 

3 y, E, 3 y, H u, 12, 13, Q, 
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of chiral fermions on the lattice M, y, lid, llJ, the analysis of the continuum limit for 

n n n 

staggered fermions 17], and the effect of quenching on the Dirac spectrum |18l. Il9|. 

For theories with a mass gap and a nonzero chiral condensate the mass dependence 
of the partition function can be obtained from a chiral Lagrangian. In particular, in the 
domain where the Compton wavelength of the Goldstone modes is much larger than the 
size of the box only the constant fields have to be included in the partition function so 
;hat its mass dependence simply follows from the mass term of the chiral Lagrangian 
2ol |. Based on this result, it was shown in [2^ that the Dirac spectrum has to satisfy 
consistency relations in the form of sum rules for inverse eigenvalues. These relations are 
not sufficient to determine the Dirac spectrum. However, it turned out that the complete 



ne c omp lete 



low-lying Dirac spectrum can be obtained from a modified chiral Lagrangian 
or equivalently, from chiral random matrix theory j^, 3]. In the Schwinger model, and 
also for QCD with Nf = 1, there are no low-lying excitations and the only contribution 
to the low-energy limit of the partition function is the vacuum energy. Starting from 
this result, it is possible to derive sum rules for the inverse Dirac eigenvalues of the 



Schwinger model |20|, |26|. However, to obtain a deeper understanding of these sum rules 
in the context of quantum field theory, it of isjnterest to derive them directly from the 
microscopic field theory. This was achieved in fZf^ for the simplest sum rule in the sector 
of zero topological charge. 

In this paper we give a microscopic derivation of the simplest sum rule in the sector of 
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arbitrary topological charge v. This requires a detailed understanding of scalar corre 



functions in the Schwin ger model which were derived before both in the plane j27[ and 



on the torus 



iH 29, SO.Isi 



ation 



It turns out that the clustering property of the scalar 



correlation functions alone is sufficient to derive the (gauge field averaged) sum rules. 
Below we will give a derivation based on this property. A third way to obtain Leutwyler- 
Smilga sum rules is to start from the bosonized form of the massive Schwinger model. In 
the microscopic derivation of sum rules we first derive a more informative quantity - the 
same sum rules for an arbitrary fixed external gauge field (which was also obtained for 
1/ = in The universal Leutwyler-Smilga sum rules then follow after averaging over 

the gauge fields. 

Analytical results obtained for the Schwinger model have been tested elaborately by 
means of Monte-Carlo simulations. The numerical value of the chiral condensate has been 
reproduced accurately 0, Q|. The index theorem was confirmed and the fermionic 
would-be zero modes were identified 0, Q. Leutwyler-Smilga sum rules have been ob- 



vyler 

0. 



tained both for zero and nonzero topological charge [10| . The distributions of the smallest 
eigenvalues of the Dirac operator for the Schwinger model 0, 3 Hi 3 agree with ana- 
ytical results obtained by means of chiral Lagrangians [2^ | 



l2fi| and random matrix theory 



A theory that is closely related to the Schwinger model are random Dirac fermions 
where the electromagnetic interaction is replaced by a random gauge field interaction. 
Starting with the work of Gade and Wegner 3^ the behavior of the spectral density 
near zero has attracted a great deal of attention in particular with applications to the 
quantum Hall effect in mind 0,0] • important difference with the Schwinger model 
is that in the condensed matter literature the physically interesting case is the quenched 
model, sometimes defined as the limit of zero fiavors. The consensus is that the density 
of states diverges for strong disorder as is the case for the quenched Schwinger model 



26 



37l |. For weak disorder, on the other hand, the spectral density vanishes. We 



will discuss the Dirac spectrum of this model and relate it to the chiral condensate and 
the simplest Leutwyler-Smilga sum rule. 

The organization of this paper is as follows. In an introductory chapter we introduce the 
Leutwyler-Smilga sum rules and the Schwinger model and discuss its main properties. A 
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combinatorial derivation of sum rules based on the clustering property of scalar correlators 
is given in section lllll Section IIVI contains a microscopic derivation of the simplest sum 
rule, and a derivation from the bosonized massive Schwinger model is given in section IVl 
It is explained in section IVTl that modified sum rules for two or more flavors can be related 
to modified clustering properties of the scalar correlators. A random matrix derivation 
of sum rules is given in section IVllI In this section we also study the fluctuations of the 
fermion determinant. In the first part of section IVlllI we discuss the effect of topology 
on the Dirac spectrum and the mechanism of the formation of a chiral condensate. In 
the second part of section IVIIII we discuss the behavior of the Dirac spectrum for random 
fermions. Concluding remarks are made in section HXl 



II. GENERAL FRAMEWORK 



A. Sum Rules 



In this paper we study the eigenvalues of the Dirac operator by means of sum rules 
for its inverse eigenvalues introduced in [2^ . The eigenvalues of the anti-hermitian Dirac 
operator are defined by 

D(f)k = ^Afc0fc. (1) 

In sector with topological charge u the Dirac operator has exactly u zero eigenvalues that 
are not paired. All other eigenvalues occurs in pairs ±Afc. In the sector of topological 
charge u the partition function is defined by 

Z'^(m) =m'' J c/A;det'[iD]e"^(^^\ (2) 

where the integral is over gauge fields in the sector of topological charge u weighted by 
the gauge field action S{Afj,). As indicated by the prime, 

det'[zD] = n (3) 

the zero eigenvalues are not included in the determinant resulting in the pre-factor m". 

At fixed 6 angle the partition function is given by 

00 

Z{m,e) = e"'^Z''{m). (4) 
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Below we will use two types of averages of an operator O denoted by single and double 
brackets. The first type is defined by (notice that Z^{m = 0) = Z{m = 0, 6)) 

{or = ^o(J^o) / DA^^O[A^^]det'me-'^^'^\ (5) 
and the second type by 



{{o)y 



J DA;0{A';)det'me-'^^r.). (6) 

To obtain sum rules for the inverse Dirac eigenvalues in a sector of topological charge 
u we expand the fermion determinant in powers of m, 

n^O ?i>0 '^n ni>0 

For the small mass expansion of the partition function we thus obtain 

^ - ™^ m . ™- (e g . ™- { E ■ ■ ■ 

= m^{det'mr[l + m^((j:^)\ + m^(( E + ' ' 

\ \n>0 / / \ \ ni>0 ^ni^n2 / / 

where we have used the definition (jH)) for the average in the sector of topological charge ly. 
Therefore, sum rules for the inverse Dirac eigenvalues follow from the sub-leading terms 
in the small mass expansion of the partition function in the sector of topological charge 
u. If 

m^(ao + a2m2 + ...), (9) 



Z^{m = 0) 
the simplest sum rule is given by 



1/ 



In j20[ general arguments where given that, for a gauge theory interacting with fermions 
according to the Dirac operator with a nonzero chiral condensate, this ratio is given by 

(11) 



4(A'/ + ^) 

For Nf = 1 this argument is particularly simple. The 6 dependence of the partition 
function, as can be observed from (j2I) and (jH), is obtained from the replacement m —>■ 
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■mexp{i6). Given that the vacuum energy is equal to —Re{mVT^), this results in the large 
volume partition function Q| 

Z{m,9) = 6"'^^''°'^. (12) 



By inverting (jlj we find that the partition function in sector u is given by 







= UmVE) = ^{mVEni + + ■■■). (13) 

This results immediately in the sum rule In this paper we will derive this sum rule 
by a microscopic calculation in the Schwinger model. 



B. The Schwinger Model 

In this section we give a brief review of the Schwinger Model Q which is massless QED 
in two dimensions. The Euclidean Lagrangian is defined by 

C=^-F^,-^[z^-g4]i^, (14) 

where the electric field strength is given by 

Foi = doA, - diAo, (15) 

and 

i ^ - g 4 = -ff,{idf,- gA^), (16) 
with Euclidean gamma matrices defined by 

7o = en, 7i = ^2, 75 = CTs. (17) 

The Lagrangian of this model has a chiral symmetry which is broken by the U{1) axial 
anomaly. As a consequence the "photon" becomes massive with mass given hj fi = g/ -^/tt. 
In this model, a local external charge is screened by massless fermions. Therefore, the 
asymptotic states contain no fermions, and the Schwinger boson is the only physical 
particle. The theory is super-renormalizable with a coupling constant g (which has the 
dimension of mass) that does not run. 
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The (massless) Schwinger model, being equivalent to a noninteracting gas of bosons, is 
exactly solvable and was solved many times, by different techni q ues (oper ator language, 
path integral, bosonization) and on different manifolds 



38, 



39 



40 



41 



42 



m. Also, 



many authors considered more specific features of the model, such as certain correlation 



functions at zero and at finite temperature |22l, ]2B, 12, 



determinant, zero modes, the index theorem and instantons 



55 



56 



57|, the multi- flavor Schwinger model jS^ 



59 



61 



44 



i 



49 



47 



50 



, the fermion 



51 



52 



53 



The vector potential can be decomposed as 

A^ = ~e^,d''<f) + d^X, (18) 

with the last term being a pure gauge. The topological charge, u, of the gauge fields is 
equal to the difference of the number of right handed and left handed zero modes of the 
Dirac operator and is therefore necessarily quantized. It is given by 



(19) 



By Stokes theorem the topological charge is determined by the asymptotic behavior of 0. 
One easily shows that the large distance asymptotics. 



ix) ~ — logx^ + const. 
2^ 



\x\ 



oo, 



(20) 



results in a topological charge equal to v. The boundary conditions (pn|l imply that the 
plane is compactified at infinity, for example by stereographic projection to a sphere with 
radius i? HHQ- 

Because of a particular property of the 2d Dirac algebra. 



7^75 = -T'f^iJ.ulu, 

the Dirac operator can be written in in the form 



(211 



(22) 



Using this representation one easily finds the following explicit expressions for the zero 
modes in the sector of topological charge z/: 

/- 



%l)p{x) 



p = ■ ■ ■ z/ — 1, 



X = xo ± ix\. 



(23) 
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Note that these zero modes are not normahzed. From the asymptotic behavior (j2U|) . it is 
easy to see that for p = 0, ■ ■ ■ , — 2 the zero modes are normahzable. However, the zero 
mode for p = u — 1 is not normahzable on the plane (its norm diverges logarithmically). 
It has to be included nevertheless because it is normahzable on the compactified plane 
as required by the index theorem The pre- factor 1 / \^2tt in the zero 

modes also follows from matching with compact manifolds |57l |. 

The unpaired zero modes explicitly break the Ua{^) symmetry of the partition function. 
In the massless limit the chiral condensate {ipip) comes from the zero modes in the sectors 
z/ = ±1. It has a particularly simple form on the plane [3,0, 57, 6^, 

S ^ {^^) = -^e\ (24) 



where 7 is the Euler constant. As a consequence of the Banks- Casher formula |66| the 
average density of the low-lying Dirac eigenvalues is given by SV^/vr. This also implies 
that if we take the thermodynamic limit before the chiral limit, we find a nonzero value 
for the condensate from the nonzero modes 



C. Effective Action 



Using the action (|T^ and the decomposition (|THjl . the partition function in the sector 
of topological charge v is given by 



= J D(j)''e-^^^'^^det[iD + m] 



n>0 



(25) 



Here we have used the chiral symmetry of the nonzero Dirac eigenvalues. The infinite 
product over the eigenvalues has to be regularized in the UV. This procedure is well 



known 







is 1^, lia |5 



and results in an anomalous contribution to the effective action. The result 



n = CdetATexp f d^x(l){x)A<p{x 



(26) 
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where C is an (infinite) normalization constant which drops out from all final results and 
J\f is the norm matrix of the (unnormalized) zero modes, 



J\fp, = J d^xijl{x)ijg{x). (27) 

The bosonic partition function in the sector of topological charge u is thus given by (to 
the lowest order in m) 



Z^m) =m'' J D(f)'' det ATe-^l"^"! , (2^ 



with the effective action F equal to 



m = l- I d'x<P{x)[A'-f,'A]<j>{x). (29) 



2 

This effective action F defines the propagator 

1 1/1 1 



In the coordinate space, 



(30) 



" ~ ^^o(/i|a;|) + const., (31) 

where KQ{fi\x\) is a modified Bessel function of the second kind, which is exponentially 
small at large distances. The large distance behavior of this propagator is therefore 
determined by the first term which is a massless propagator 

^(a;) = -T^logx^ + const. . (32) 



The arbitrary constant in the equation ()31|) can only be obtained after a suitable 
infrared regularization of the theory. For exam ple, thi s constant can be determined by 
matching with the result on compact manifolds 



iple, this 
0,0- 



D. Partition Function in the Sector of Topological Charge 



To evaluate the partition function in the sector of topological charge u (which we take 
positive for notational convenience; of course, the partition function does not depend on 



11 



the sign of the topological charge), we first rewrite the norm matrix. Using the explicit 
expressions for the zero modes we find 

detAr= ^ / d2xi---d2a;,l^sg(a7r)xr^°^xi^°^---xr^'^"'^x;:(""'^e-'^^ 

^ ■ (JIT 

The sums over permutations a and vr can be rewritten as the product of two determinants, 

2 



u-l 



1 Xi ■ ■ ■ Xi 



I/y ... nr'^ 



-1- / d'x, ■ ■■d'x.e-'^^U'''^^^^ n - 



X 



(34) 



where a well-known property of Vandermonde determinant was used. Using this expres- 
sion in (j21), we find for the partition function to the lowest order in m. 



(35) 



Since the effective action is Gaussian the path integral is simply given by 

ZO(m = 0) 7 ^ ■ V ; 

Using the explicit expression for Q on the plane ()3ip we observe that the factor 

ni>7 is canceled by the asymptotic behavior of the Greens function. Thus, 



in the limit of large volume 



m 



Z^{m = 0) I/! V 27r 
The chiral condensate can be expressed as 

1 



1 , Z^im) + Z-^im) 
lim 



Z^{m = Qi)Vm^o m 
Since Z^^[m) = Z^{m), the chiral condensate is given by 



le232g(0) 
77 



(37) 



(3^ 



(39) 



which can be used to eliminate ^(0) from the partition function. Naively, eliminating ^(0) 
in favor of the condensate is equivalent to fixing the constant in ()3ip. Notice that only the 
product of the square of the normalization constant of the zero modes and exp(2(7^^(0)) 
is fixed by this condition. 
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III. COMBINATORIAL DERIVATION OF SUM RULES 

In the sector of topological charge v we find for the small mass expansion of the 
partition function in a finite but large volume V 

^0, = E ^ / ^'^1 ■ ■ ■c^'^n(^^(a;i) ■ ■ (40) 

Z^\m = 0) nl J 

The average is over field configurations with topological charge u weighted by the gauge 
field action and the determinant of the Dirac operator in the space of nonzero modes. 
For one fiavor, the terms in the expansion with n < v vanish, since the number of zero 
modes in the correlator is not sufficient to compensate the zeros of fermion determinant. 
For n = u we have 

For massless quarks, this expectation value is called a minimal correlator and saturated 



by zero modes [^Ol, |4^ 1^ . If we use the explicit expressions for the zero modes and take 
into account their (anti-)symmetrization, we recover (jH^ and what followed. However, 
since we are only interested in the leading large-volume behavior, there is a simpler way 
- to use the clustering property: for all | 

{^rM^i) ■ ■■^rM^u))" ^ (|)', (42) 



so that, in agreement with (jHTjl 



Let us consider sub-leading orders in m in (jlUj) . For n > u, the non- vanishing contribu- 
tions have n — V contractions from the nonzero mode part of the propagator. Since the 
massless propagator connects only states with the same chirality the expectation value 
(i/jiplxi) ■ ■ ■il)il){xyj^i)Y = 0. The first non-vanishing sub-leading term in is 

— — — / (fxi ■ ■ ■ d^Xn(fxd^y{ijij{xi) ■ ■ ■ ijij{xn)iJij{x)ijij{y)y. (44) 
[u + 2)1 J 

This correlator is sometimes called a non-minimal correlator and was calculated in 



68|. However, for our purpose we will only need the large distance limit of this correlator. 
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which again follows from the clustering property of the correlators. By decomposing the 
field into left handed and right handed components we find 

(f xi. . .(f XydP x(f y {ipTp^xi) ■ ■ ■ 'ip'ip{xy)'il)4'{x)'il)'ilj{y)) y (45) 
= (z/ + 2) / d^xi ■ ■ ■ d^Xyd^xd^y{'4}LiiR{xi) ■ ■ ■ ^LijR{xy)i}Li)R{x)i)Ri)L{y)Y ■ 



The combinatorial factor z/ + 2 arises because the fermion bilinear with opposite chiralities 
can be at each of the v + 2 points. Using the clustering property for the expectation value 
we obtain for the large volume limit of (j4Uj) 



+ 2X2) =7^(2) • («) 



(z/ + 2)!^ '^2^ (z/ + 

The mass dependence of the partition function to this order is thus given by 



Z0| 

Comparing this result with (jH)), we arrive at the sum rule 



1 1 y2T/2 



which also incorporates the case of negative v. In j20[ this result was obtained from the 
mass dependence of the partition function at nonzero vacuum angle. In our derivation we 
only used the clustering property of expectation values. Therefore our derivation is also 
valid for QCD with one massless fiavor. For the Schwinger model a derivation of the sum 

n 

rule for zero topological charge was given in |26| . 

The term of the order nf^'^\ I = 2, 3,... is not much more complicated. Again 1/ of the 
fermion bilinears are saturated by zero modes. The remaining ones are contracted by the 
nonzero mode Green's function, and therefore we have to choose / of the z/ + 2/ fermion 
bilinears with opposite chirality j69|, 



J d^xi ■ ■ ■d^Xy+2i{i'i'{,xi) ■ ■ ■'iIj'iIj{xu+2i)Y (49) 

d^Xi ■ ■ ■ d'^Xy+2l{^L^R{Xl) ■ ■ ■ ^pL^R{Xu)tljL^R{x,y+i ■ ■ ■ ^Ri)L{Xu+2l)Y ■ 




Again using the clustering property we obtain for the large volume limit of the m"^"^^ 
order term in (jlUj) . 

(z/ + 2n! /El^x"+2/ ^v+2i , , 



(z/ + 2/)! (z/ + /)!Z! V 2 ^ (z/ + /)!/!V 2 
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The term of order m^' in the expansion of ((7|) is given by 

1 



\2 . . . \2 

ni>0 ^rii 



(51) 



Gathering all the pieces, we obtain the final result 



20 







I 


2' n (|z/ 


1+0! 



2i 



(52) 



This sum rule is valid for any theory with the above clustering property. In particular, it 
is valid for QCD with one flavor 20|. 



IV. MICROSCOPIC DERIVATION OF SUM RULES IN THE SCHWINGER 

MODEL 



A microscopic derivation of the sum rule (j4Hj) in the sector of topological charge u = 
was given in |22i]. In this section we generalize this derivation to arbitrary nonzero 
topological charge. The main idea of the derivation may be summarized in the following 
identity 

E 4 = -Tr[G'^'], (53) 

where is the Green's function of the Dirac operator with external field = —enpdp(j). 
The equality (|53|) follows from the spectral representation of the Green's function 

g>.,).i:^^(£m:m. ,54) 

and the fact that excited Dirac states are normalized to 1 (the trace is both over Dirac 
indexes and spatial coordinates). Therefore, we need the Green's function G^{x^y) for an 
arbitrary background field in the sector of topological charge v (which is taken z/ > for 
convenience). Due to the index theorem, the Dirac equation has exactly v right-handed 
zero modes for a gene ric field configuration in this class. The Green's function G'^{x,y) 
is thus defined by |70| 

I p.^Cix, y) = 5{x - y) - itP'i^^ y)^ (55) 
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where 



p=i ! d?z'iljl{z)i)p{zy 
is the projector on the subspace of zero modes, ipp{x), p = 1, ■ ■ ■ , z/, and 

75^ = ^(1+75). 



(56) 



(57) 



In other words, Gu{x,y) is the Green's function in the space of the nonzero modes. The 
exphcit solution of (fH^ is given by: 

G''{x,y) = (1 -75+P'^)G^(l -75+n, (58) 
where, because i 1/)^ = e^'^'^^i ^e^'^^, 

G'ix^y) = e-3^(")^^G'o(x - y)e-3'^^y^^' . (59) 



(60) 



The free two-dimensional Dirac propagator Go{x,y) is given by 

Go{x,y) = ^ ^f'''y'\ 
2tt [x — yY 

In the representation ipi{x) = e~^'^^^^ and {ipk, A; = 2, ■ ■ ■ , i^} is an orthogonal set perpen- 
dicular to ip, one can easily show that (|KH|l is a solution of Using that •y^G'^'j^ = 
we find the sum rule 

1 



-Tr 



-Tr 



;i-(i + 75)n[Gt 
;i + 75)(i-n[G'i' 



(61) 



In the second equality we have used that in the term proportional to the identity, both 
helicities give the same result. Using the explicit representation for the Green's function 
and carrying out the trace over the 7 matrices, we find 

1- = - I d^xcfyd^iSix -y)- P'^(x, i/))tr [(1 + 75)e-^^^(2/)75^^(^^ ^)g-230{.)75^^^^^ ^)g-90{x)75 
(fx(fy(fz{5{x -y)- p'^(a;, y))e-9<^(!^)-9<^(^)+29<^W 



X 



(2vr)2 

(x - z)^(z - y)f, - ie^^jy - z)f,{z - x)^ 
{x — zY{z — yY 

d^xd^ydh{6{x - y) - P''(x, |/))e-^'^(^)-^'^(^)+^f'^(^) 

1 1 



X 



(27r)2 

x-yf - 2i€^^{y - z)^{z - x 



{x — zY{z — yY 



(x — zY (y — zY 



(62) 
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In the last equality we have used that 



(x - z)^{z - y)^ = h{x - yf - (x - zf - {z - yf). 



(63) 



The last two terms in the square brackets are independent of either y or x and vanish 
after integration as can be seen from the representation ()56|) 



d'yP''{x,y)e 



(64) 



Using that the integral over the 5-function also vanishes, we find the sum rule 
11" 



A, 



(2vr)^ 



d xd yd zP^{x, y)e 



g4>{y)~g4>{x)+2g<t>{z) ~ V) ~ 2^6^^ (y — 2:)^(2: — x)^ 

(x — zy{z — yY 



(65) 



The projector formula (j^UI) is written in terms of orthogonal modes. However, since 
the zero modes ()23p are not orthogonal, we require a projector formula that is valid for 
an arbitrary set of non-orthogonal modes. This is easily achieved by using a Lagrange 
interpolation such that the projector is equal unity for each of the zero modes. The 
required form is given by 



P\x,y) 



^sgn( 



(66) 



where the normalization factor is the familiar determinant of the norm matrix, given by 
fl34|l . The projector can be expressed in terms of Slater determinants as 

/ d'^X2 ■ ■ ■ d'^Xyx{x, X2,---, Xy)x\y, a;2, ■ ■ ■ , Xy) 



P''{x,y) = u^- 



(67) 



J d'^xi ■ ■ ■c?2x^|x(xi, ■ ■ ■ ,x,,)|2 

Here, x{xi, ■ ■ ■ , Xu) is the Slater determinant of fermionic zero modes. In this representa- 
tion, it is particularly transparent that the projector is not only independent of the choice 
of basis, but (unlike det A/"), is also independent of the normalization of the zero modes, 
as of course it should be. Using the explicit representation of the fermionic zero modes 
we find the projector 

Jdx2 ■ ■■dxyWk=2ix+ - xt){y- - X, )A({x+, ■ ■ ■ ,x+})A({x2 , ■ ■ ■ , x;})e-^(^(-)+<^(J^))-25(^(-^)- 



(^.)) 



V 



Jdxi--- dXyA{{xt, ■ ■ ■ , x+})A({xr, ■ ■ ■ , x-})e-'^9(<f>{^^)+-+'f>M) 



(68) 
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where A is a Vandermonde determinant. Renaming x by Xi and y by x^^i, the sum rule 
in the sector of topological charge u can be written as 



TTT = 7 ^ — TTTT Tr"; 77 f d? Z [ dXi ■ ■ ■ dXy+iC 

At^o^n (27r)-+>-l)!detAry J 



2g{,i,(xi)+-+,i,{x^+i))+2gct>{z) 



X 



k=2 



{xi - zy{z - Xu+iy 

(69) 

Only the symmetric part of the expression in the second line of this equation contributes 
to the integrand. Since the expression is already symmetric in X2, Xj^, we only need to 
symmetrize with respect to Xi and x^+i. Then we get, 

V 

^5 n (^^ - ^^fc - ) A({a;^, ■ ■ ■ a;+})A({x2 , ■ ■ ■ x;}) 

k=2 

(xi - zf'{z - 

+ 1) ^5 Ui^piXp - XT) Ui^g{x+ - xl){xp - zf{z - Xg)^' 

(70) 

The term with p = q vanishes and can be included in the sum. Using that 

p Y\k^p{Xp Xk) 

we can replace the numerator 

{xp - Xqf/2 - ie^^{xp - z)^{z - Xq)y (72) 

in dZOl) by 

{Xp Xq^ /2 {Xp z) /2 {Xq z) /2 i^^^viXp z) ^(^Z Xq)i/ 

= -(x^ - z+)ix- Z-). (73) 



After simplifying the numerator we find for (f7n|) 
2A({x+ •••a;++J)A({xr,---a;;+i}) 



E 



1 Ui^piXp -Xi )Ui^q{x+ -xl){x- -z-){x+ -z+)' 

(74) 
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The sum in this equation factorizes into the product of two sums. Each sum can be 
expressed determinant 



E 



A({xi,---x^^i}) 



i/(xr 



1 



- v~\ 



X 



- u-1 
u+l 



• • X 



}) 



(75) 



rip (-^p ^ , 

The determinant is easily evaluated by multiplying each row by a factor (x^ Exactly 
the same simplification can be made for the + variables. The sum rule therefore simplifies 
to 

1 2 



(27r)^+2(z/+ 1)! det 



■2g{0(xi)+-+(/.(x^+l))+2g0{^) 



A({a;+ 



a; 



})A(K 



■ a; 



}) 



n,^il(^+-x+)(z--x-) 



(76) 



^iipdp4> in the sector of 



This sum rule is valid for a arbitrary external gauge field = 
topological charge u. The Leutwyler-Smilga sum rules (fTn|) are obtained after averaging 
over the fields (p{x) with the effective action det A/'e^^''^^, where r[0] is defined in (jSHI)- 
Notice that effective action in the nontrivial topological sectors, due to the presence of 
det A/", is not only non-Gaussian in (p but also non-local in (p. However, the factor det A/" 
cancels in the average of the sum rule. The resulting path integral is Gaussian in and 
the average over the fields is given by a lowest order cumulant expansion. This results 
in 



E \2 



A{{xf, ■ ■ - x 



X 



})A({xi,---x^+i}) 



X' 



The second line of (f77j) is canceled by the asymptotic behavior of the exponentiated 
Green's functions in the first line of ()77|) (see ()32|)). The large volume limit of the sum 
rule is therefore given by 



(77) 



E \2 



(271)^^+2(2/+ 1)! 



yv+2^2(u+2)g^G(Q) 



19 



where we used (jSH))- Now, using (|^. and the definition (jH)), we finally obtain for the 
Leutwyler-Smilga sum rule. 



E 



m=0 \A„7^0 



A? 



2 u 



SV. (79) 



(27r)-+2(z/+ 1)! 

The result ()79p is in agreement with the sum rule obtained by Leutwyler and Smilga 
and with the result obtained from Random Matrix Theory 24 1 . 



V. BOSONIZATION AND SUM RULES 



In this section we derive sum rules for the inverse Dirac eigenvalues starting from a 
bosonic description of the massive Schwinger Model. In the particle physics literature, 
abelian bosonization goes back to the work of Lowenstein and Swieca |2|, who essen- 
tially bosonized massless Schwinger Model, and Coleman I7ll and Mandelstam [2^. The 
bosonized form of the massive Schwinger Model was studied in |73l . I74| . Most earlier 
works used an operator language. The path integral approach to bosonization, which we 
will use below, was developed later, in many works, in particular in [3 Hi- 
lt is well-known [tsI . ItJ ] that the action of the massive Schwinger model can be written 

as 



S[^] 



(fx 



I 2 ^f 



^0) 



with the constant c = fie'^'^ /{27i) = The value of this constant is such that the chiral 
condensate is equal to its known value (j2H)- At small m this is a theory of a weakly self- 
coupled massive scalar field. The physics of this model is discussed in detail in 3, Iz^ | . 
In the mean field limit, where = 0, we obtain ()12p from the action (jHU)). This directly 
leads to the Leutwyler-Smilga sum rules. 

The first step in the microscopic derivation of sum rules is to eliminate the mass term 
of the ip field by 



(81) 
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After partial integration of {dipY and shifting (/? by we obtain the partition function 

Z{m,e) = j D^pD(t)e-^^«^'^''^\ (82) 



with effective action given by 

-1 



<S'efr[0,¥'] 



d X 



g2g0+i 2^if+ie _|_ g-2g</>-i 2^tf- 



(83) 



The integral is only over modes with A0 not identically equal zero. The partition 
function in the sector of topological charge v is given by the Fourier transform of the 
partition function, 

Jo 27r 

The integral over 9 is easily performed by expanding the integrand in powers of e*^, 



X 



9=1 



k,l=0 
I 

p=i 



(85) 



It is clear that only the terms with k — contribute to the integral. However, we can 
also perform the integral after shifting 



In the effective action this results in the extra term 

ie^ j d^xAcf). 

The integral over 9 therefore vanishes unless 

i.e. the topological charge of the fields is equal to u. 
The path integral over ip is Gaussian, and is given by 



(86) 



^7) 



(88) 



Dip 



lie 

p=i 



i 2^/TT(f{yp) 



v+l 

He 

q=l 



-i 2^/n<p{xq) 



n'^gtinUiyp - 



(89) 



with 
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So[^] = \j (fx{d,ipf. (90) 



Thus we end up with 

Z^{m) = J D^e-^«[-l / D0e-^[^1 £ ^^^^^^ J d% ■ ■ ■ d\d'x, ■ ■ ■ d'x.^, 



1=0 + 



The path integral over can be performed by a leading order cumulant expansion. The 

Vandermonde determinants are again canceled by the large distance asymptotic behavior 

of the Green's functions. In the large volume limit we thus obtain 

Mm) ^ ^^'^ 
Zo(m = 0) fr"^ (z/ + /)!/!V 2 ^ ' ^ ' 

which is equal to the result obtained from the clustering assumption The sum rules are 

the same as derived in section ED 

n n 

From this result one easily derives that Z{m,9) = exp(mV^S cos 0) (see also [20,|22(). 

We end this section with a conjecture for sum rules in a fixed external 0— field. By now 
we realize that the field introduced in ()81|) is indeed the same as in previous sections. 
We can extract the sum rules from comparison with the partition function in eigenvalue 
representation. To that end, we rewrite the partition function as 

J det 

OO 1 

e-r[<^W^detAr^5:m2' ^ (94) 

1=0 "iT^O '^ni ' ' ' ^ni 

where the formulas (|7j) and (j^Tj) have been used. We remind that deiM'^ is determinant 
of the norm matrix (j34jl . We conjecture that the integrands of (j9T| and (j94|l are equal: 



\2 



<jK-^t -^3) ^2g\<t>{yi)+-+cl>(yi)-<)>{x^) 0(x,+i)] 



For / = 1 this conjecture was shown to be valid in in section HVl 
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VI. CLUSTERING FOR TWO OF MORE FLAVORS IN CHPT 

In multi-flavor QCD, cliiral symmetry is broken spontaneously, with the appearance 
of massless Goldstone bosons. This results in long range correlations that will modify 
the clustering property of scalar correlators and will affect the sum rules for the inverse 
Dirac eigenvalues. Starting from a chiral Lagrangian we will propose a modified clustering 
relation. 

Let us first show that factorization leads to an incorrect result. To leading order in 
the quark masses the partition function in the sector of topological charge v is given by 

^XXnXZ''' = ''"'1.!^^'^" / n ■ • • «'(4") • • ■ • ■ ■ **(4''")>". 

(96) 



If we assume clustering as in the case of one flavor we obtain 



vNf 



resulting in the partition function 



(9^ 



ZO(mi = 0) [z/!]^/ 

However, instead of , the correct pre-factor j2^ is given by Ilfc^o ^ ^'/('^ + ^)'- 

The reason for the incorrect pre-factor in ()98p is that we did not account for the 
spontaneous symmetry breaking of SUA{Nf), which results in massless Goldstone bosons 
and long range correlations. As a consequence, the clustering property does not hold 
directly even for the correct ground state (i.e. at fixed 9). It is possible, as we will show 
below, to write down a modified clustering property. 

The clustering property follows from the large distance behavior of the correlators. 
This is determined by the zero momentum term in the chiral Lagrangian, i.e. by 

with normalization constant A/" such that the partition function is equal to unity for 
M = 0. For a diagonal mass matrix we obtain the expansion 

^Low I 



(100) 
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which should be compared to the large volume limit of the QCD partition function given 
by 

no Nf 



fci,---,fcjVy=0 1=1 



(101) 



We thus obtain the clustering relation 



Nt 



Nf 



(102) 



1=1 JU^SU(Nf) 

If we isolate the phase of the diagonal matrix element Uu this can be rewritten more 
suggestively as 

Nf Nf 

{I[i^i^i)'')o= [ , ,'iUVL\Mm%iNf)-i,^ (103) 

JV(iSV{Nf) \.^l"f> ?i 

with {4'4'){e/Nf)-^i the chiral condensate for one flavor and vacuum angle equal to {9/Nf) — 
^i. Notice that only the diagonal matrix elements of U enter in this expression. Since 
\Uu\ < 1 we could interpret the factors \Uii\ as the projection of the one flavor condensates 
onto certain directions in group space. 

In general, the integral in p02|) cannot be evaluated analytically. However, the situ- 
ation simplifies enormously for the minimal correlators in the sector of fixed topological 
charge u. In this case we have rii = n2 = ■ ■ ■ = n^^ = u and the integral over 6 is only 
nonzero if no factors Ul occur. This results in 



Nf 

1=1 



2 







vNf , Nf 

duxiui 

uesu{Nf) 



(104) 



This integral can be evaluated analytically (see appendix A) 



n 

i=0 



(105) 



The leading mass dependence of the partition function is obtained by substituting the 
result ()104|) into equation (j^ . 



,mNf 





V 




_ 2 




2 



n 



1=0 



I + V) 



(106) 



This explains the suppression of the scalar correlator for N ^ > 2. 
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VII. RELATIONS WITH RANDOM MATRIX THEORY 

The chiral random matrix model in the sector of topological charge v and Nt flavors 

nn 

with equal quark mass m is defined by j2j, |25|| 

Zl^{m) = J dWdet^f{D + m)e-"^'^^^\ (107) 

where 

D=\ \ L {m 




and ly is a complex n x (n + i/) matrix. The eigenvalues of this model are distributed 
according to a semicircle. The largest eigenvalue is equal to 1 /S. We use this eigenvalue to 
normalize the partition function such that it becomes dimensionless. If we also normalize 
the partition function to unity for Nf = we obtain the normalization constant 

/ \ n{n+u) 

C(i) = e"^M-) (109) 



We also have included the factor e^^f which eliminates the constant vacuum energy 
for Nf 7^ 0. We can evaluate the partition function ()107|) in the limit m — ^ in two 
different ways. First, using an eigenvalue representation of W , and second, using a a- 
model representation of the large n limit of the partition function. Since the cr-model 
is exactly the partition function (jHS)), we only discuss the first approach in the next 
subsection. 

A. Eigenvalue Representation 

To obtain an eigenvalue representation of the partition function we use the polar de- 
composition 

W = UAV-\ (110) 

where U G U{n + z/)/f/"(l) x f/(z/) and V G U{n). The Jacobian of this transformation 
is given by 

n 

J = A({A^})2"n(A.)'''^^ (Ill) 

k=l 
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The eigenvalues of A are positive or zero, and the integration is over an ordered sequence 
of eigenvalues which will be accounted for by a factor l/n\ below. To lowest non- vanishing 
order in m, the eigenvalue representation of the partition function is thus given by 



/ n 2rfA.Af ^^^)^^A^({A^})e-«^^S.^s (112) 



k=l 

where we have introduced a second constant 



^(2) ^ g-n2_^^ vol {U (n + z/))vol {U {n)) 
vol([/(z/))vor(f/(l)) ' 



(113) 



and the volume of U{n) is given by 



(27r)' 



n(n+l)/2 



Yo\{U{n))= ^ ■ (114) 

llfc=o 

We rewrite the partition function in terms of the integration variables 

Xk = 2nT,'^\l (115) 

This results in 

1 r 

Z'N,{m) = C^^C^^m^^f^^ / n dx,x:^^A\{x,})e-^^l\ (116) 

k=i 

where 

(^(3) _ ^2i2)-["("+^)+"^/]7]-~"("+'^)e"'^-''. (117) 

The integral over the Xk was calculated by Forrester ^\ . Our final result for the partition 
function is thus given by 

ra-l 

Z];^,(m) = m''^SY.-'^nNffj{Z)fjm^n(u+Nf+n) ^-j^^ + Nf+ (118) 

j=0 

Collecting all factors we find the partition function 

(m) = ^-A^/s-2"^/n-"^^e"^^ TT ('^ + ^/ (ng) 

Using the Stirling formula, we find in the thermodynamic limit 

Z;^,(m) = m-">S-'""r '''' . (120) 
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For Nf = 1 we find the ratio 



7U=1 

'Nt=l 



7U=0 



(171) 



m{n + l). (121) 



According to the condensate is given by 

1 



hm hm 

y— >oo m-*o mV 



(122) 



In random matrix theory the volume is identified as = 2n such that ()121|) results in a 
condensate equal to unity. The reason for this incorrect result is that we have evaluated 
the ratio of two products with the same number of nonzero eigenvalues. To obtain a result 
with the correct dimensionality we have to regularize the determinants by only including 
eigenvalues below a given energy 0]. For = 1 we have on average one eigenvalue less 
with absolute value below a given energy than for z/ = 0. This results in an overall factor 
of 1/pA, where p is the number of eigenvalues below the cut-off energy and A is the 
spacing between the eigenvalues. Instead of using a fixed energy cut-off, we calculate the 
ratio of the partition functions for a fixed total number of eigenvalues meaning that in 
the sector of topological charge v we make the replacement n n — v/2 m p2()|l . We 
finally find the ratio 

— [mnhj ^ II 



(123) 



in agreement with general arguments given in |2C)| . This result also follows from the fact 
that the large limit of the random matrix model is given by the nonlinear sigma model 
[24j . As was already found in ,24j , the result p23j) confirms that the chiral anomaly 
is contained in random matrix theory. 

Since the finite volume partition function for Nf = 1 is a Bessel function, it is tempting 
to identify its zeros with the average position of the eigenvalues. In the random matrix 
model it can be easily checked whether this identification is justified. Numerically, it 
turns out that the ratio obtained from the average eigenvalues is a factor 2/tx smaller 
than the result given in p21|l . A second question we have asked is whether the average 
determinant is mainly given by the product of the average eigenvalues or whether it is 
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due to the fluctuations of tlie eigenvalues. It turns out tliat tlie product of tlie average 
eigenvalues (with the determinant included in the weight) is an order of magnitude larger 
than the average determinant. 



VIII. DISCUSSION 

A. Spectral Duality 

The chiral condensate can be obtained in two ways. First, by taking the chiral limit 
after the thermodynamic limit, 

S(^) = lim lim :^/t^^ ), (124) 



m^o v^oo V \ iXk + m 
and, second, by inverting the two limits. 



= lim lim ^ ). (125) 



V^oom^O V \iXk + fri/ 
In the first case, the condensate arises as a consequence of spontaneous symmetry break- 
ing, and in the second case it is due to the anomaly or instantons. The two condensates 
are not necessarily equal. For example, if we restrict the partition function to gauge field 
configurations with zero topological charge, the second definition gives zero. By spectral 



duality we mean that the two condensates are equal for the 6 - vacuum [79|. In the sec- 
ond case, when the contribution to the condensate comes from the sector with topological 
charge equal to one, the eigenvalues for = 1 are shifted such that the condensate is equal 
to the one obtained from the first definition. The first definition implies the Banks-Casher 
formula 6^ 

S« = ^. (126) 
This implies that the average spacing of the eigenvalues is given by 

A = = —. (127) 

p(o) ^ ' 

If we assume that the position of the eigenvalue A„ is given by nA, we obtain the sum 
rule 

y. = -sv^, (128) 
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which explains the functional dependence of the sum rule. 

As we will show next, the large v limit of the sum rule can be obtained from the large 
V limit of the average position of the eigenvalues. In the sector of topological charge v 
the "average" position of the eigenvalues can be expressed in terms of the zeros, j^^i^, of 
the Bessel functions 



A. = ±1^. (129) 



For large k at fixed v one one finds 



Afc~j:;,A~(fc + ^-5)A with 5=1. (130) 

Here, A = l/p(0) is the average spacing of the eigenvalues. This formula does not 
accurately give the position of the eigenvalues for large v and finite k. In this case one 
can use a uniform asymptotic expansion of the zeros given by 

jV,fc ~ yzk, (131) 



with Z]^ implicitly defined by [Ak — l)7r = Au{^ z1 — 1 — arccos(l/zfc)). This results in the 
correct large v limit of the sum rule 

E^-SV^/ dk'=^! dz^^ = ^. (132) 



Let us now calculate the condensate according to the second definition 

^(2) ^ 1 (nl(^Afc + m))^=i ^ 1 (nfc(^Afc + m))^=_i 
V {Y{k{i^k + ni))^=o V" (nfc(^Afc + m))^=o ■ 



(134) 



Using the expression ()T3|) for the partition function, this can be rewritten as 



= hm - 

*o V 



m- 



Ii{mVIl) /m ^ /_i(myS)/m 



loimVJ:) /o(m\/S) 

The ratio in ()133p can be interpreted as the ratio of the expectation values of the fermion 
determinant in the topological sector u = 1 and z/ = 0. 

We will calculate the averages in p33j) by replacing the eigenvalues according to p29j) 

"'-'OV {Uk{i^k + m))^=o ™-^oV^nfc(^J^=o,fc/(S\/)+m) 

^ 1 U'k U.=i,fc A/tt Uk U^Lo.fc A/tt U'k A/tt 
V U'k ^JZi,kN^ Uk U.=o,fc A/vr Uk U."Lo,fc A/^ ■ 
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The first two ratios are finite, but the third ratio has to be regularized which we will do 
by a function regularization (see Appendix B). For m = we find the result (see 

1 m U.^Li,,A/7r 1 T\l-6) S«r2(l-5) 



The ratios Y['k'^ju,k/Y['k'^ju=i k "^^^ be easily evaluated numerically. We find 



(136) 



Uk iju=o,k 

Hk ^j"=0,fe 



Wk i3u=i,k 

Y\k '^iu=i,k 



TC 



r2(1.25)/2V2. 



(137) 



Therefore, 



(138) 



Another way of evaluating the ratio 



is from the limit 



lim 



1 m=-n^jZl,k^/^ 
nUk=-n^jT=0,kNT^ 



lim — 

n — ^oo fi 



T-rn-1 
llfe=0 



A' i 



1/2 



(139) 



Using the infinite product representation of the F-function we find that this limit is given 



by F2(|)/F2(|). We thus have 



1 1 nLoU^Li,feA/7r 
— iim 



(140) 



A™nnLo^&,fcA/7r' 

so that the definition ()139j) of the infinite product would have resulted in an incorrect 
value for the chiral condensate. This explains why in the random matrix calculation of 
previous section the incorrect chiral condensate is obtained from 1)1211) . The fact that the 
two results differ by the eigenvalue at the cut-off is in agreement with the interpretation 
that C-function regularization includes the product of the eigenvalues up to a fixed energy. 



B. Random Gauge Field 



In the condensed matter literature Dirac fermions in random gauge fields have received 

a considerable amount of attention. Among the different types of models that have been 

I I 

considered, we mention the random magnetic field problem iSQl, the random flux model 

n n n 

|35| . the random mass model j3^ and the random gauge field model j3^. In this section 
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we only consider the latter model which is the quenched Schwinger model with the gauge 
field action is replaced by 

1 „o 



2^2 M M 



The factor g'^ is introduced to have the same normalization as in 
we can make a Hodge decomposition of the gauge field [sil. 

Let us consider the action of an instanton configuration given by 

1 



(141) 
Also in this case 

(142) 
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\og{x + p 



(143) 



In the Schwinger model the action of this configuration is given by 



TT 



(144) 



which is infrared finite. Also the anomalous contribution to the effective action in the 
Schwinger model is infrared finite, 



rf2^^<^a^^=l(l0gp2 + l). 



(145) 



For the random gauge model, on the other hand, we obtain the action 



o / ^ *^^// — n / 

2 J\x\<R ^ 2 cr^ j\x\<r 



TT 



2cr2 



log( 



2 



P 



R^ + p"^ 



(146) 



which diverges for R oo which implies that instantons are suppressed in the thermo- 
dynamic limit. Therefore the chiral condensate of the massless theory is zero, and by 
spectral duality, the eigenvalue density of the Dirac operator vanishes independent of the 
topological charge. For cr^ = vr the suppression factor is 1/V, which exactly as in the 
Schwinger model for Nf = 1, results in a finite spectral density for E 0. Notice that 
the difference between the anomalous contribution in effective action of the Schwinger 
model and the random gauge action is only in the boundary term, which diverges for 
R ^ oo. 
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Let us evaluate the spectral density of a random gauge model with Nj massless flavors 
and n flavors with mass m. Because of the above remark, we restrict ourselves to the 
trivial topological charge sector, where the boundary term that gives the diverging result 
in ()146|) is absent. Using that the pure gauge term decouples from the partition function, 
the effective action of this model is given by 

J ^ Z ZTT Q J 

The massless fermion fields result in an overall constant which can be ignored. The action 
of the massive fermion fields will be replaced by their bosonized form given by 

-j^ n n 

S'f(0,(/9„) = - ^{d^Lfaf -cmY^ co?,{2ipa^fn -2ig(j)). (148) 
^ 0=1 0=1 

If we define 

r(0) = -^0(4^' + /^'(^/ + - (149) 

the total bosonic action is given by 

^RG = j d^x[T{<P) + Sf{<P, ^o)]- (150) 



After shifting the fields by 



we obtain the action density 



e , , 

- 77^, (151) 



1 n 1 2 n n 

^ E + o + ^'^Nf)i^^(|>f -cmY, cos(2(^,v^ -9)- ifid,<l) ^ d^^Va- (152) 

^ a=0 ^ ^ OL=\ a=\ 

The final action density is obtained after performing the Gaussian integral over 0, 

n n 1 " ^ 

oC EE d^,ipe,d^i{>i3 + - ^{d^ipaf - cm E cos{2ipaVT^ - 0), (153) 

^ a=l /3=1 ^ a=l a=l 

where 

C = (154) 
Considering the cosine term as a vertex, the propagator of the fields is given by 



Gc.M = -,i^o.p--^). (155) 
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Using this result one can easily derive the dependence of the chiral condensate on the 
infrared cutoff which is taken to be equal to the size of the box. Following the derivation 
of js^] and using the lattice spacing a as ultraviolet cutoff we find for ajh that 

Go.,{\x\ - 0) = - — ^) log J (156) 

The chiral condensate is given by 

E = ii9„logZ~(e^^v^^i) 
n V 

In the replica limit, n — >^ 0, we obtain the condensate 

1 



(158) 



According to the Banks-Casher the smallest nonzero eigenvalue can be estimated as 

A„,„ ~ ^ ~ (159) 
On the other hand, if the eigenvalue density is given by 

p(A) ~ VX", (160) 
the smallest nonzero eigenvalue follows from 

r 

Jo 

Therefore, 



p{X)dX ^ 1. (161] 



so that 



a = - — i. (163) 

1 + C 

This result has also been derived from a more sophisticated renormalization group analysis 
83[ . For the action (|153|) we obtain 



82 



^ ^^ + /iV(iV,-l) _ l + a^(iV,-l)/7r 
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In the quenched case {Nf = 0), this result agrees with the result obtained in l34ll. and 
in the limit a ^ oo we recover the result for the Nf flavor Schwinger model [2 (J . IgsI i . 
For Nf = 1 the spectral density is qualitatively different from the Schwinger model. A 
diverging spectral density for Nf = has been observed in recent lattice simulations . 
If the spectral density vanishes, it still possible to derive sum rules for the inverse Dirac 



eigenvalues 
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It could be that such sum rules are easier obtained for lattice 
'ermions than for lattice simulations of the Schwinger model 



simulation of random Dirac 
with two massless flavors [18 1. 

A special case where the sum rules for the inverse Dirac eigenvalues can be evaluated is 
Nf = and n = 1. Since the sum rule follows from the coefficient of in the expansion 
of the partition function in powers of m, the spectral density corresponding to this sum 
rule is the one of the Schwinger model with one massless flavor. 

For n = 1 and Nf = the partition function after the transformation ipi —>■ /i0 is given 

by 



(165) 



where 
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Let us evaluate the simplest sum rule in the sector of zero topological charge. To this end 
we expand the partition function to second order in m and integrate over 6. We obtain 

dl\ogZ^=°{m) = J ^2x^6'"'''''''^"'"^-'"'''''°^°'°^ (167) 
where the random gauge fleld Green's function is given by 

^^^(x, y) = log(x - y f. (168) 

The volume dependence of the sum rule is thus given by 

9^1ogZ^=°(m) ~ V^"^'/^\ (169) 
If the eigenvalue density scales as 

p{E) ~ VE", (170) 
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the volume dependence of the sum rule is given by 



e4~^^, (171) 

k ^k 



so that the exponent a is equal to 



Indeed, this result agrees with (jl64p for Nj = 1. 



IX. CONCLUSIONS 



In conclusion, we have analyzed Leutwyler-Smilga sum rules in the Schwinger model 
and found that they are in agreement with the universal result from chiral perturbation 
theory and random matrix theory. Instead of relying on general symmetry arguments, we 
have performed a microscopic calculation of the sum of the inverse square Dirac eigenval- 
ues in the sector of arbitrary topological charge z/, generalizing a result by Smilga obtained 
in the sector of zero topological charge. 

We have shown that validity of the sum rules follows from the clustering property 
of the scalar correlation functions. This argument also applies to QCD with one flavor. 
For QCD with several flavors, the naive clustering argument fails due to the presence of 
Goldstone bosons. A modified clustering property is obtained from the chiral Lagrangian 
that corresponds to the low energy limit of QCD partition function. 

The dependence of the sum rules on the topological charge v is consistent with a shift 
of the Dirac eigenvalues by v/2 times the average level spacing. Such shift of the Dirac 
spectrum exactly results in the chiral condensate for massless quarks in the sector of 
topological charge u = 1. However, obtaining this result, requires a regularization of 
the fermion determinant with a fixed energy cutoff such as for example the (^-function 
regularization. 

In the microscopic derivation, the sum rules where obtained from more general sum 
rules valid for a fixed external gauge field with topological charge u. The most general 
sum rules in this class were conjectured based a simplified derivation of Leutwyler-Smilga 
sum rules starting from the bosonized Schwinger model. It would be interesting to probe 
such sum rules directly in lattice simulations. 
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The Schwinger model is part of a larger class of models known under the name of 
random Dirac fermions which have been investigated in the context of the quantum Hall 
effect. The key difference between the two models is that instantons are suppressed by 
the random gauge field action. Therefore, the chiral condensate of the massless theory 
vanishes in the domain below a critical value of the disorder. In these models we expect 
modified sum rules that would be a useful probe for the scaling behavior of the density 
of states near zero. 

Acknowledgements. We would like to thank Andreas Ludwig, Andreas Smilga and 
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X. APPENDIX A 

In this appendix we calculate the integral given in ()1()5|1 . We start from the result 

W{J)= I dU exp{Tr{JU)) 

JsU{Nf) 

We are interested in the integral 

'/ d Y f d 



dUTF Tl" ■ ■ ■ TF 

SUiN,)^''^^^ ^^^^^ 



\dJiiJ \dJNj,Nf 



WiJ) 



(174) 
j=o 



On the other hand, we can choose J = diag[Jii, J22, ■ • • , JNfNf], so that 

det J = J11J22 ■ ■ ■ ^AT^AT^. (175) 

Using this result, we arrive at Nf independent differentiations over diagonal matrix ele- 
ments of J. The only term in the sum that contributes is the one with i = u resulting 
in 

Nf 2!---(A/-l)! 



dUU^,,U^,---U^^^ 

SU{Nf) ^ ^ 



u\ ■ ■ ■ {u + Nf - 1 

Nf-l ., 

,1 ^^f 



n (176) 



i.o (' + !') 
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XI. APPENDIX B 

In this appendix we compute the determinant of the Dirac operator in ^-function 
regularization with eigenvalues given by the asymptotic values of the zeros of Bessel 
functions, 



j^k-k+'^-S with 5^^. 



(177) 



In a C function regularization the determinant of the Dirac operator in the sector of 
topological charge v is given by 



log^'(^AJ:y7^ + 



d_ 
ds 



s=0 



E 



+ E 



^i{k + ^-5 + im/AyA' t^i(k + ^-S- im/AyA 



(178) 



The function in between the brackets is known as the Hurwitz zeta function, 

oo 2 

For convenience we only calculate the product for m — 0. For u — 1 we find 



[C(s, - - 5)A-] 

s=0 ^ 

= -2C'(0,^-5) + 2C(0,^-5)logA. 



Combining this with the result for i/ = 0, 



d 



Us,\-b)A- 



= -2C'(0,l-5) + 2C(0,l-5)logA, 



(179) 



:i80) 



(181) 



and using that 



we find the ratio. 



C(0,a) = \-a, C'(0,a) = log r(-i log(27r). 



i nl^fjgii., ^ 1 r^(i-^) ^ r^(i-^) 



(182) 
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